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Linear response methods are applied to identify the increase in electron-phonon coupling in el- 
emental yttrium that is responsible for its high superconducting critical temperature T c , which 
reaches nearly 20 K at 115 GPa. While the evolution of the band structure and density of states is 
smooth and seemingly modest, there is strong increase in the 4d content of the occupied conduction 
states under pressure. We find that the transverse mode near the L point of the fee Brillouin zone, 
already soft at ambient pressure, becomes unstable (in harmonic approximation) at a relative vol- 
ume V/V = 0.60 (P w 42 GPa). The coupling to transverse branches is relatively strong at all high 
symmetry zone boundary points X, K, and L. Coupling to the longitudinal branches is not as strong, 
but extends over more regions of the Brillouin zone and involves higher frequencies. Evaluation of 
the electron-phonon spectral function a 2 F(ui) shows a very strong increase with pressure of coupling 
in the 2-7 meV range, with a steady increase also in the 7-20 meV range. These results demonstrates 
strong electron-phonon coupling in this system that can account for the observed range of T c . 



I. INTRODUCTION 

The remarkable discover in 2001 of MgB2 with 
superconducting critical temperature T C =40K, and 
the fact that the simple free-electron-like metal 
lithium^ also has T c in the 14-20K under 30-50 
GPa pressure, has greatly increased efforts in seeking 
higher T c in elements and simple compounds. Cur- 
rently 29 elements are known to be superconducting 
at ordinary pressure and 23 other elements supercon- 
duct only under pressure^. Among elements there 
is a clear trend for those with small atomic number 
Z to have higher values of T c , although much vari- 
ation exists. For example, under pressure^ Li, B, 
P, S, Ca, and V all have T c in the range 11-20 K. 
Hydrogen2*i2iii, the lightest element, is predicted to 
superconduct at much higher temperature at pres- 
sures where it becomes metallic. 

While light elements tend to have higher T c 
among elemental superconductors, Hamlin et al. 12 
recently reported that Y (Z=39) superconducts at 
T C =17K under 89 GPa pressure and 19.5 K at 115 
GPa, with the trend suggesting higher T c at higher 
pressure. This result illustrates that heavier ele- 
ments should not be neglected; note that La (Z=57) 
has T c up to 13 K under pressurei 13 i 14 The super- 
conductivity of La has been interpreted in terms of 
the rapidly increasing density of states of 4/ bands 
near Fermi level with increasing pressure, causing 
phonon softening and resulting stronger coupling un- 
der pressurei 15 ' 16 Such a scenario would not apply 
to Y, since there are no / bands on the horizon 
there. No full calculations of the phonon spectrum 
and electron-phonon coupling have been carried out 
for either Y or La to date. 



La and Y are two of the few elemental transition 
metals to have T c aboveZ^ 10 K, and the case of Y is 
unusually compelling, since its value of T c is at least 
as high that of Li, qualifying it as having the highest 
T c of any elemental superconductor. Moreover, the 
reduced volume v = V/Vn=0A2 corresponds to the 
value of T c w 20K in Y 12] (115 GPa) and also to 
the report of T c 20K in (strained) Li 0] above 50 
GPaii2*i& For our study of Y reported here, it is first 
necessary to understand its phase diagram. Under 
pressure, it follows a structure sequencei2*2°i through 
close-packed phases that is typical of rare earth met- 
als: hcp^Sm-type— >dhcp^dfcc (dfee is distorted 
fee, with trigonal symmetry). The transitions occur 
around 12 GPa, 25 GPa, and 30-35 GPa. Supercon- 
ductivity was first found— in Y by Wittig in the 11- 
17 GPa range (1.2-2.8 K), in what is now known to 
be the Sm-type structure. From 33 GPa to 90 GPa 
T c increases smoothly (in fact T c increases linearly 
with decrease in v over the entire 35-90 GPa ranged) 
suggesting that Y remains in the fee phase, perhaps 
with the distortion in the dfee phase vanishing (the 
tendency is for the c/a ratio in these structures to 
approach ideal at high pressure^). Calculations^ 
predict it adopts the bec structure at extremely high 
pressure(>280 GPa), but this is far beyond our in- 
terest here. 

In this paper we report electronic structure and 
electron-phonon coupling calculations of Y for re- 
duced volumes in the range 0.6< v < 1 (pressures up 
to 42 GPa). Our results indeed show strong electron- 
phonon coupling and phonon softening with increas- 
ing pressure. A lattice instability (in the harmonic 
approximation used in linear response calculations) 
is encountered at v—Q.Q and persists to higher pres- 



sures. The instability arises from the vanishing of 
the restoring force for transverse displacements for 
Q||< 111 > near the zone boundary, correspond- 
ing to sliding of neighboring close-packed layers of 
atoms. It is only the stacking sequence of these lay- 
ers that distinguishes the various structures in the 
pressure sequence of structures (see above) that is 
observed in rare earth metals. Near-vanishing of 
the restoring force for sliding of these layers is con- 
sistent with several stacking sequences being quasi- 
degenerate, as the structural changes under pressure 
suggest. 

This paper is organized as follows. In Sec. II 
structural details are given, and the calculational 
methods are described. Results for the electronic 
structure and its evolution with pressure are pro- 
vided in Sec. III. The background for understand- 
ing the electron-phonon coupling calculations is pro- 
vided in Sec. IV, and corresponding results are pre- 
sented and analyzed in Sec. V. The implications are 
summarized in Sec. VI. 



II. STRUCTURE AND CALCULATION 
DETAILS 

Yttrium crystallizes in the hep structure at am- 
bient pressure with space group P63/mmc (#194) 
and lattice constants a=3.647 A and c=5.731 A^i 

Since the observed structures are all close packed 
(or small variations from) and above 35 GPa Y is 
essentially fee, we reduce the calculational task by 
using the fee structure throughout our calculations. 
The space group is Fm3m (#225), with the equiv- 
alent ambient pressure lattice constant a=5.092 A. 
We do note however that results for electron-phonon 
strength can be sensitive to the crystal symmetry, 
both through the density of states and through the 
nesting function that is described below. 

We use the full potential local orbital (FPLO) 
coda^i to study the electronic structure, and ap- 
ply the full-potential linear-mumn-tin-orbital (LM- 
TART) code 2 ^ to calculate the phonon fre- 
quencies and the electron-phonon coupling spec- 
tral function a 2 F. For FPLO, a k mesh of 
36 3 and the Perdew-Wang (PW92)2& exchange- 
correlation potential are used. The basis set is 
Is2s2p3s3p3d::(4s4p)/5s5p4d+. For LMTART, a 
k mesh of 48 3 and GGA96 approximation 27 for 
exchange-correlation potential are used. For the 
electron-phonon coupling calculations we used a 
phonon Q mesh of 16 3 , which has 145 Q points in 
the irreducible Brillouin zone. 
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FIG. 1: (Color online) Plot of the total DOS and pro- 
jected Ad DOS per atom of fee Y with different volumes. 
Both the total and the Ad density of states at Fermi level 
decrease with reduction in volume. 



III. ELECTRONIC STRUCTURE UNDER 
PRESSURE 

Many studies suggest that the general character of 
an elemental rare earth metal is influenced strongly 
by the occupation number of the d electrons, which 
changes under pressure. Our calculations show that 
the Ad occupation number of trivalent Y increases 
from 1.75 at ambient pressure, to a little above 2 
at V=0.7Vo and then finally close to 3 at V=0.3Vo 
(which is extreme pressure) . Such an increase can be 
seen from the projected density of states (PDOS) of 
4c? states (Fig. ^) at different volumes. From Fig. ^ 
broadening of the density of states with reduction in 
volume can be seen, but is not a drastic effect. The 
main occupied Ad PDOS widens from 2 eV to 3 eV 
with reduction of the volume to v—0.5. The value of 
of the density of states at the Fermi level (taken as 
the zero of energy) N(0) decreases irregularly with 
volume reduction; the values are given in Table [I] 

The pressure evolution of the band structure is in- 
dicated in Fig. where the Ad character at v=1.00 
(black) and t>=0.50 (gray) is emphasized. First, the 
relative positions of the Fermi level crossings change 
smoothly, indicating there is little change in the 
Fermi surface topology. This slow change is also 
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FIG. 2: (Color online) Plot along high symmetry direc- 
tions of the bands of Y at V/V o =1.00 and at V/V o =0.50. 
The "fattening" of the bands is proportional to the 
amount of Y 4d character. Note that the 4d character 
goes substantially in the occupied bands under pressure 
(the lighter shading), although there is relatively little 
change in the Fermi surface band crossings. 



seen in Fermi surface plots, of which we show one 
(below). Second, the overall band widths change 
moderately, as was noted above in the discussion of 
the density of states. The change in position of 4d 
character is more substantial, however. Ad bands at 
X lying at -1 eV and -2 eV at ambient volume are 
lowered to -3 eV and -4 eV at v=0.50. Lowering of 
4<i character bands at K and W is also substantial. 
Thus Y is showing the same trends as seen in al- 
kali metals. For Cs and related alkalies and alkaline 
earths under pressure, 6s character diminishes as 5c? 
character grows strongly with pressured In Li, 2s 
character at the Fermi surface evolves to strong 2p 
mixture^ at the volume where T c goes above 10 K. 

The Fermi surface of Y at ambient pressure (hep) 
has been of interest for some time, from the pioneer- 
ing calculation of Loucks^ to the recent measure- 
ments and calculations of Crowe et al~ However, 
the unusual Fermi surface in the hep structure (hav- 
ing a single strong nesting feature) is nothing like 
that in the fee phase we are addressing, which is 
unusual in its own way. At u=1.00 the fee Fermi 
surface is a large 'belly' connected by wide necks 
along < 111 > directions as in Cu, but in addition 
there are tubes ('wormholes') connecting a belly to 
a neighboring zone's belly through each of the 24 W 
points. The belly encloses holes rather than elec- 
trons as in Cu; that is, the electrons are confined to 
a complex multiply-connected web enclosing much 




FIG. 3: (Color online) Surface plot of the Fermi surface 
of fee Y at a volume corresponding to ambient pressure. 
The surface is shaded according to the Fermi velocity. 
The surface is isomorphic to that of Cu, except for the 
tubes through the W point vertices that connect Fermi 
surfaces in neighboring Brillouin zones. The evolution 
with pressure is described in the text. 



of the surface of the Brillouin zone. 

As the volume is reduced, the wormholes slowly 
grow in diameter until in the range 0.5< v <0.6, 
they merge in certain places with the necks along 
the < 111 > directions, and the change in topology 
leaves closed surfaces around the K points as well as 
a different complex multiply-connected sheet. The 
point we make is that, at all volumes, the Fermi 
surface is very complex geometrically. There is little 
hope of identifying important "nesting" wavevectors 
short of an extensive calculation. Even for the sim- 
ple Fermi surface of fee Li, unexpected nesting vec- 
tors were located 34 in three high symmetry planes 
of the zone. The rest of the zone in Li still remains 
unexplored. 



IV. BACKGROUND: ELECTRON-PHONON 
COUPLING 

The electron-phonon spectral function a 2 F(u>) 
can be expressed in terms of phonon properties [and 
N(0)] in the form2fi 



a 2 F(uj) = 



J2^6(w-w Qv ) (1) 



2ttA^(0) — lu q „ 



where 



F(w) 



Qi> 



(2) 



is the density of phonon states, and N(0) is the single 
spin Fermi surface density of states. The phonon 
linewidth 7q„ is given by 

7Q „ = \M£l Qtk \ 2 6(su)S(e k+Q ) (3) 



where M^_q is the electron-phonon matrix element; 
v is the branch index. Sums over Q or k are conven- 
tionally normalized (divided by the number of unit 
cells in the normalization volume). 

The quantities thus defined enable one to identify 
the contribution to A from each mode, the "mode 
A" , as 



A 



0J r 



N(0) 



k 



\ 2 S(e k )S(e k+Q ) (4) 



With this definition A is the average over the zone, 
and sum (not average) over the N v — 3 branches of 
all of the Xq u values. The electron-phonon coupling- 
strength A then is given by 



A = 



7rJV(0) 



E 



— = E 

Qi/ 



(5) 



The critical temperature T c can be obtained by 
using the Allen-Dynes modification of the McMillan 
formulaj2iii2£ which depends on the logarithmic, first, 
and second frequency moments u>i ogi ojj =< u> >, 
and uj2 =< to 2 > 1 / 2 , as well as A and the Coulomb 
pseudopotential fi* . These averages are weighted ac- 
cording to the normalized coupling 'shape function' 
2a 2 F (oj) / (Xoj) . They are often, and will be for Y es- 
pecially under pressure, much different from simple 
averages over the spectrum F(o>). 

Note that Aq^, or 7q„, incorporates a phase space 
factor, the 'nesting function— £(Q) describing the 
phase space that is available for electron-hole scat- 
tering across the Fermi surface(EF=0), 



^Q) = 4E <5 ( £ k)^k +Q )cx I 

1_ w C 



c |vk x « k +Ql 

(6) 

Here C is the line of intersection of an undisplaced 
Fermi surface and one displaced by Q, and Vf. is the 
electron velocity at k. These equations presume the 
adiabatic limit, in which the phonon frequencies are 
small compared to any electronic energy scale. This 
limit applies to elemental Y. 




FIG. 4: (Color online) Plot of the calculated phonon 
spectrum along high symmetry directions (F-X, T-K, 
T-L) of fee Y with different volumes. The longitudi- 
nal mode phonons increases with the distance from V 
points along all the three directions. Along F-X direc- 
tion (left panel), the doubly degenerate traverse mode 
slightly softens near X point, while along F-K direction 
(left panel, only the T2 mode sightly softens near K 
point. Along F-L direction (right panel), the already 
soft doubly degenerate transverse mode soften further 
near the L point with decreasing volume. At V = O.6V0, 
the frequency at L becomes negative, indicating lattice 
instability. 



V. RESULTS AND ANALYSIS 
A. Behavior of Phonons 

The calculated phonon branches are shown along 
the high symmetry lines, from v—0.90 down to 
v—0.60, in Fig. 0] The longitudinal modes behave 
normally, increasing monotonically in frequency by 
~ 30% in this range. The transverse modes along 
< 100 > and < 110 > show little change; the dou- 
bly degenerate transverse mode at X softens by 20%, 
reflecting some unusual coupling. Along < 1 10 >, 
Ti and T2 denote modes polarized in the x — y plane, 
and along the z axis, respectively. 

The interesting behavior occurs for the (doubly 
degenerate) transverse branch along < 111 >. It is 
quite soft already at t>=0.9 (7 meV, only 25% of the 
longitudinal branch), softer than the corresponding 
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FIG. 5: (Color online) Plot of the calculated linewidths 
of fee Y for varying volumes. The linewidths of the trans- 
verse modes at the X point increases from 1.3 to 5.5 
as volume decreases from V=0.9Vo to V=0.6Vo. The 
linewidths of the T2 along < 110 > modes show the 
same increase. The linewidths along the < 111 > direc- 
tion have been multiplied by four for clarity. 



mode in hep Y at ambient pressured With decreas- 
ing volume it softens monotonically, and becomes 
unstable between v=0.65 and v=0.6Q. It should 
not be surprising that the transverse mode at the L 
point is soft in a rare earth metal. The sequence of 
structural transitions noted in the introduction (typ- 
ically hep— >Sm-type— >dhcp— >dfcc— >fcc for trivalent 
elements) involves only different stacking of hexago- 
nal layers of atoms along the cubic (111) direction. 
So although these various periodic stackings may 
have similar energies, the soft (becoming unstable) 
transverse mode at L indicates also that the bar- 
rier against sliding of these planes of atoms is very 
small. At v=0.60 (see Fig. the largest instability 
is not at L itself but one-quarter of the distance back 
toward T. At w=0.65 there are surely already anhar- 
monic corrections to the lattice dynamics and cou- 
pling from the short wavelength transverse branches. 
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B. Linewidths 

The linewidths 7q„, one indicator of the mode- 
specific contribution to T c , are shown in Fig. [S] 
To understand renormalization, one should recog- 
nize that in lattice dynamical theory it is lu 2 , and 
not u> itself, that arises naturally. At v=0.90, u> 2 
for the T modes is only 1/16 of the value for the 
longitudinal mode at the L point. A given amount 
of coupling will affect the soft modes much more 
strongly than it does the hard modes. 

For the < 110 > direction, the strong peak in jq v 
for the T2 z polarized) mode at the zone bound- 
ary point K (5.7 meV) is reflected in the dip in 
this mode at K that can be seen in Fig. At 
u=0.60 the linewidth is 1/3 of the frequency. The 
coupling to the Ti mode along this line is negligible. 
Note that it is the Ti mode that is strongly coupled 
in Li and is the first phonon to become unstable. 
A peak in the linewidth of the L modes correlates 
with a depression of the frequency along this line. 
Along < 001 > the T modes again acquire large 
linewidths near the zone boundary under pressure. 
This electron-phonon coupling is correlated with the 
dip in the T frequency in the same region. 

The coupling along the < 111 > direction is not 
so large, either for T or for L branches (note, they 
have been multiplied by a factor of four in Fig. 
|S| The coupling is strongest at the zone boundary, 
and coupled with the softness already at w=0.90, 
the additional coupling causes an instability when 
the volume is reduced to v=0.60 (P = 42 GPa). 
This seems to represent an example where a rather 
modest amount of coupling has a potentially catas- 
trophic result: instability of the crystal. Evidently 
Y is stabilized in the fee structure by anharmonic ef- 
fects, coupled with the fact that being already close- 
packed there may be no simple structural phase that 
is lower in energy. 

C. Coupling Strength 

It is intuitively clear that strong coupling to ex- 
tremely low frequency modes is not as productive in 
producing high T c as coupling to higher frequency 
modes. This relationship was clarified by Bergmann 
and Rainer^ who calculated the functional deriva- 
tive ST c /5a 2 F(ui). They found that coupling at fre- 
quencies less than Q = 2-7rT c has little impact on T c 
(although coupling is never harmful). Since we are 
thinking in terms of Y's maximum observed T c w 
20 K, this means that coupling below Q = 10 meV 
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FIG. 6: (Color online) Plot of the product Xq v luq v of 
fee Y for different volumes, along the high symmetry 
directions. Note that the longitudinal (L) values along 
< 111 > have been multiplied by four for clarity. In 
addition, values corresponding to unstable modes near L 
have been set to zero. Differences in this product reflect 
differences in matrix elements; see text. 



becomes ineffective. 

The product XqvWqv oc Jqu/wqu gives a some- 
what different indication of the relative coupling 
strength^! than does either \q v or 7q„. It is also, up 
to an overall constant, just the nesting function de- 
fined earlier, with electron- phonon matrix elements 
included within the sum. Since the nesting function 
is a reflection of the phase space for scattering, it is 
independent of the polarization of the mode, hence 
differences between the three branches are due solely 
to the matrix elements. 

This product \q v ioq v is shown in Fig. The 
weight in the transverse modes is concentrated near 
the zone boundary, with the region being broader 
around L than at X or K and growing in width with 
pressure. The T\ branch along < 110 >, which 
is polarized along [110], shows essentially no cou- 
pling. The weight in this product for the longitudi- 
nal modes is peaked inside the zone boundary along 
the < 001 > and < 110 > directions with a mean 
value of 7-8 meV, and is weaker along the < 111 > 
direction. 



TABLE I: For each volume v studied, the columns give 
the experimental pressure™^ GPa), the Fermi level den- 
sity of states N(0) (states/Ry spin), and calculated val- 
ues of the mean frequency uj 1 —< lu > (meV), the 
logarithmic moment uji og and second moment uj-z =< 
uj 2 > 1//2 (all in meV), the value of A, the product Xui 2 
(meV 2 ), and T c (K). For T c the value of the Coulomb 
pseudopotential was taken as /x*=0.15. 
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7.9 


6.9 


9.5 


12.1 


2.80 


410 


(16.9) 



D. a 2 F(uj) 

The results for a 2 F are displayed in Fig. The 
longitudinal peak in the 20-35 range hardens nor- 
mally with little change in coupling strength. The 
7-20 meV range of transverse modes at v—0.90 in- 
creases in width to 2-24 meV at w=0.60, and the 
strength increases monotonically and strongly. The 
strong peak in a 2 (uj) = a 2 F \ui) j 'F '(w) , shown in the 
bottom panel of Fig. reflects the very soft modes 
that have been driven into the 2-5 meV range, and 
the fact that they are very relatively strongly cou- 
pled. The substantial increase in coupling, by a fac- 
tor of ^2.5, in the range 7-25 meV is important for 
T c , as noted in the next subsection. 



E. Estimates of T c 

This background helps in understanding the 
trends displayed in Table|U where T c from the Allen- 
Dynes equation 3 ^ (choosing the standard value of 
/i*=0.15) and the contributing materials constants 
are displayed. The calculated values of A increases 
strongly, by a factor of 3.7 in the volume range 
we have studied. Between i>=0.65 and v=0.60 (the 
unstable modes are removed from consideration) 
A increases 30% but T c increases by only 2.5 de- 
grees. The cause becomes clear in looking at the 
frequency moments. These moments are weighted 
by a 2 F(uj)/uj. a 2 (uj) itself become strongly peaked 
at low frequency under pressure, and it is further 
weighted by The frequency moments set the 

scale for T c (uji og in particular) and they decrease 
strongly with decreasing volume. In particular, 
ojiog decreases by 45% over the volume range we 
have studied, reflecting its strong sensitivity to soft 
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to A, but is also very effective in lowering the tem- 
perature scale (ujiog). For a 2 (u>) shapes such as 
we find for Y, the quantities A and uji og which go 
into the Allen-Dynes equation for T c do not provide 
a very physical picture of the change in T c . For 
this reason we provide also in Table Q] the product 
Xui 2 = N(0) < I 2 > /M (< I 2 > is the conventional 
Fermi surface average of square of the electron-ion 
matrix element and M is the atomic mass). For the 
volumes 0.60< v <0.80 in the table, the ratio of 
Xu 2 /T c is nearly constant at 23±1.5 (in the units of 
the table), illustrating the strong cancellation of the 
increase of A with the decrease in frequency moments 
in producing the resulting T c . 



VI. SUMMARY 
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FIG. 7: (Color online) Top panel: Plot of a 2 F(uj) versus 
lu As volume decreases, a 2 F(uj) increases and gradually 
transfers to low frequency. Bottom panel: the frequency- 
resolved coupling strength a 2 (uj) for each of the volumes 
studied. The evolution with increased pressure is dom- 
inated by strongly enhanced coupling at very low fre- 
quency (2-5 meV). 



modes. 

The increase in T c probably owes more to the in- 
crease in coupling in the 10-25 meV range (see a 2 (uj) 
plot in Fig. El a factor of roughly 2.5) than to the 
more spectacular looking peak at very low frequency. 
Put another way, the very low frequency peak in a 2 F 
looks impressive and certainly contributes strongly 



In this paper we have presented the evolution of 
elemental Y over a range of volumes ranging from 
low pressure to 40+ GPa pressure (V/V — 0.60). 
Lattice instabilities that emerge near this pressure 
(and persist to higher pressures) make calculations 
for smaller volumes/higher pressures unrealistic. For 
simplicity in observing trends the structure has been 
kept cubic (fee); however, the observed phases are 
also close-packed so it was expected that this re- 
striction may still allow us to obtain the fundamen- 
tal behavior underlying the unexpectedly high T c 
in Y. On the other hand, the Fermi surface geom- 
etry varies strongly with crystal structure, and the 
nesting function £(Q) and perhaps also the matrix 
elements may have some sensitivity to the type of 
long-range periodicity. 

In addition to the band structure, Fermi surface, 
and electronic density of states, we have also pre- 
sented the phonon dispersion curves and linewidths 
along the high symmetry directions, and also have 
presented a 2 F(cu) and the resulting value of T c . The 
results show indeed that Y under pressure becomes 
a strongly coupled electron-phonon system, readily 
accounting for value of T c in the range that is ob- 
served. 

In spite of having used a relatively dense mesh of 
Q points for the phonons, it seems clear that this 
Brillouin zone integral is still not well converged. 
Evaluation of £(Q) on a very fine Q mesh in three 
planes for fee Li, which has a very simple Fermi sur- 
face, has shown^i that this nesting factor contains 
(thickened) surfaces of fine structure with high inten- 
sity. The convergence of this zone integral (and for 
example the resulting a 2 F function) has rarely been 
tested carefully in full linear response evaluations of 



phonons; such a test could be very computationally 
intensive. Nevertheless, the general finding of strong 
coupling is clear. 

Very recently it has been found that isovalent Sc is 
superconducting at 8.1 K under 74 GPa pressured 
Note that if the lattice were harmonic and the only 
difference between Sc and Y were the masses (which 
differ by a factor of two), T c = 20 K for Y would 
translate to T c = 28 K for Sc. (For an element 
with a harmonic lattice, A is independent of mass.) 
The corresponding argument for (again isovalent) La 
gives T c = 16 K. La has T c = 13 K at 15 GPa, and 
has not been studied beyond^ 45 GPa. 

Another comparison may be instructive. Dynes 
and Rowell obtained and analyzed tunneling data 42 
on Pb-Bi alloys where A is well into the strong cou- 
pling region, becoming larger than two as is the case 
for Y under pressure in Table [I] The Pbo.65Bio.35 
alloy has A=2.13, < u? >=22.6 meV 2 . We can 
compare directly with the w=0.65 case in Table [I] 
which has A=2.15, < uj 2 >=159 meV 2 . The prod- 
uct M < uj 1 > for Y is three times as large as for 
the heavy alloy. Since the A's are equal, the value 
of N(0) < I 2 > (equal to AAf < uj 2 >) is also three 
times as large as in the alloy. The values of T c are 
14.4 K (Y) and 9 K (alloy) [somewhat different val- 
ues of //* were used.] The values of uJi og differ by 
less than a factor of two, due to the low-frequency 
coupling in a 2 (to) in Y that brings that frequency 
down, and that is why the values of T c also differ by 
less than a factor of two. 

While this study is in some sense a success, in 
that it has become clear that strong electron-phonon 
coupling can account for the remarkable supercon- 
ductivity of Y under pressure, there remains a seri- 
ous shortcoming, one that is beyond the simple lack 
of numerical convergence that would pin down pre- 
cisely A, T c , etc. What is lacking is even a rudimen- 
tary physical picture for what distinguishes Y and 
Li (T c around 20 K under pressure) from other ele- 
mental metals which show low, or vanishingly small, 
values of T c . 
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The rigid muffin-tin approximation (RMTA) of 
Gaspari and Gyorffy^ which approximates the 
phonon-induced change in potential and uses an 
isotropic idealization for the band structure to de- 
rive a simple result, seemed fairly realistic for the 
electronic contribution (the Hopfield nf) for transi- 
tion metal elements and intermetallics42*^i On top 
of these idealizations, there is an additional uncer- 
tainty in < u 2 > that must be guessed to obtain A 
and T c . One would not 'guess' the values of the fre- 
quency moments that we have obtained for Y under 
pressure. 

In addition, the RMTA expression does not distin- 
guish between the very different matrix elements for 
the various branches, giving only a polarization and 
Fermi surface average. Nevertheless, it gave a very 
useful understanding of trends^ in electron-phonon 
coupling in elemental transition metals and in some 
intermetallic compounds. While the linear response 
evaluation of the phonon spectrum and the result- 
ing coupling seems to work well, this more detailed 
approach has not yet provided - even for elemental 
superconductors - the physical picture and simple 
trends that would enable us to claim that we have 
a clear understanding of strong coupling supercon- 
ductivity. 
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